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We investigate the possibility of exciting self-oscillation in a perpendicular ferromagnet by the 
spin Hall effect on the basis of a nonlinear analysis of the Landau-Lifshitz-Gilbert (LLG) equation. 
In the self-oscillation state, the energy supplied by the spin torque during a precession on a constant 
energy curve should equal the dissipation due to damping. Also, the current to balance the spin 
torque and the damping torque in the self-oscillation state should be larger than the critical current 
to destabilize the initial state. We find that these conditions in the spin Hall system are not satisfied 
by deriving analytical solutions of the energy supplied by the spin transfer effect and the dissipation 
due to the damping from the nonlinear LLG equation. This indicates that the self-oscillation of a 
perpendicular ferromagnet cannot be excited solely by the spin Hall torque. 
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I. INTRODUCTION 

Nonlinear dynamics such as fast switching and self¬ 
oscillation (limit cycle) has been a fascinating topic in 
physics!*!. Magnetization dynamics excited by the spin 
transfer effect!^ in a nanostructured ferromagnet! - — 
provide fundamentally important examples of such non¬ 
linear dynamics. The magnetization switching was first 
observed in Co/Cu metallic multilayer in 2000!. Three 
years later, self-oscillation was reported in a similar 
system!. In these early experiments on the spin transfer 
effect, linear analysis was used to estimate, for exam¬ 
ple, the critical current destabilizing the magnetization 
in equilibriu m 13 1 14 . However, recently it became clear 
that nonlinear analysis is necessary to quantitatively an¬ 
alyze the magnetization dynamic s 2 ' 15- — . For example, 
current density to excite self-oscillation can be evaluated 
by solving a nonlinear vector equation called the Landau- 
Lifshitz-Gilbert (LLG) equatio n 23 ' 24 . 

Originally, the spin transfer effect was studied by ap¬ 
plying an electric current directly to a ferromagnetic mul¬ 
tilayer. Recently, however, an alternative method em¬ 
ploying the spin Hall effect has been used to observe the 
spin transfer effect— -4 °. The spin-orbit interaction in a 
nonmagnetic heavy metal scatters the spin-up and spin- 
down electrons to the opposite directions, producing a 
pure spin current flowing in the direction perpendicular 
to an applied current. The pure spin current excites the 
spin torque, called spin Hall torque, on a magnetization 
in a ferromagnet attached to a nonmagnet. The direc¬ 
tion of the spin Hall torque is geometrically determined^!, 
and its magnitude shows a different angular dependence 
than the spin torque in the ferromagnetic multilayer—. 
Therefore, it is fundamentally unclear whether the phys¬ 
ical phenomena observed in the multilayer— - !! can be re¬ 
produced in the spin Hall system, and thus, new phys¬ 
ical analysis is necessary. The magnetization switching 
of both in-plane magnetized and perpendicularly mag¬ 
netized ferromagnets by spin Hall torque was recently 
reported!! -31 ' 36 !’ 7 . Accordingly, it might be reasonable 


to expect reports on self-oscillation by spin Hall torque. 
However, whereas self-oscillation has been observed in 
the in-plane magnetized system—, it has not been re¬ 
ported yet in the perpendicularly magnetized system. 

The purpose of this paper is to investigate the possibil¬ 
ity of exciting self-oscillation by spin Hall torque based on 
a nonlinear analysis of the LLG equation. We argue that 
two physical conditions should be satisfied to excite self¬ 
oscillation. The first condition is that the energy that the 
spin torque supplies during a precession on a constant en¬ 
ergy curve should equal the dissipation due to damping. 
The second condition is that the current to balance the 
spin torque and the damping torque in the self-oscillation 
state should be larger than the critical current to desta¬ 
bilize the initial state. This is because the magnetization 
initially stays at the minimum energy state, whereas the 
self-oscillation corresponds to a higher energy state. We 
derive exact solutions of the energy supplied by the spin 
transfer effect and the dissipation due to damping in the 
spin Hall system by solving the nonlinear LLG equation, 
and find that these conditions are not satisfied. Thus, the 
self-oscillation of a perpendicular ferromagnet cannot be 
excited solely by the spin Hall torque. 

The paper is organized as follows. The physical condi¬ 
tions to excite a self-oscillation is summarized in Sec. El 
These conditions are applied to the spin Hall system in 
sec. uni Section El is devoted to the conclusions. 

II. PHYSICAL CONDITIONS TO EXCITE 
SELF-OSCILLATION 

Let us first summarize the physical conditions neces¬ 
sary to excite self-oscillation. The magnetization dynam¬ 
ics are described by the LLG equation 

dm TT rr / \ dm 

-= -7m x H 7-ffs m x (p x m) + am x -, ( 1 ) 

at at 

where m and p are the unit vectors pointing in the 
directions of the magnetization and the spin polariza- 
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tion of the spin current, respectively. The gyromag- 
netic ratio and the Gilbert damping constant are de¬ 
noted as 7 and a, respectively. The magnetic field H 
relates to the energy density of the ferromagnet E via 
H = -dE/d(M m), where Ad is the saturation magneti¬ 
zation. The strength of the spin torque, H s , is propor¬ 
tional to the current density j. Since the LLG equation 
conserves the norm of the magnetization, the magnetiza¬ 
tion dynamics can be described as a trajectory on a unit 
sphere. The energy density E shows constant energy 
curves on this sphere. For example, when the system has 
uniaxial anisotropy, the constant energy curves are lati¬ 
tude lines. The self-oscillation is a steady precession state 
on a constant energy curve excited by the field torque, 
the first term on the right-hand side of Eq. (|TJ). This 
means that the second and third terms of Eq. © , aver¬ 
aged over the constant energy curve, cancel each other. 
In other words, the energy supplied by the spin trans¬ 
fer effect during the precession on the constant energy 
curve equals the dissipation due to the damping. This 
condition can be expressed as^^ 

/ rl F 1 

dt— = W S + W a = 0 , ( 2 ) 

where the energy supplied by the spin transfer effect and 
the dissipation due to the damping during the precession 
on the constant energy curve of E are given 

W S (E) = 1 M dtH s [p H - (m p) (m H)], (3) 


W a (E) = —ajM 


dt 


H“ 


(m-H) 2 


( 4 ) 


The time integral is over a precession period on a con¬ 
stant energy curve. We emphasize that Eqs. m and © 
are functions of the energy density E. We denote the 
minimum and maximum values of E as E m - m and E max , 
respectively. When the energy density also has saddle 
points E s addle, E'max in the following discussion can be 
replaced by E sa( To excite the self-oscillation, there 
should be a certain value of the electric current density 
that satisfies Eq. © for E min < E < E max in a set of 
real numbers. Therefore, Eq. © can be rewritten as 


3 j g K, W s + W a =0. (5) 


We denote the current satisfying the first condition, Eq. 
©. or equivalently Eq. ©, as j(E). 

Another condition necessary to excite self-oscillation 
relates to the fact that the magnetization initially stays at 
the minimum energy state. To excite any kind of magne¬ 
tization dynamics, the spin torque should destabilize the 
initial state, which means that a current density larger 
than the critical current density, j c = j(E min ), should be 
injected. Then, the condition 

3(E) > j (-Emin) 7 ( 6 ) 



FIG. 1: (a) Schematic view of system. The current density 

j flows in the nonmagnet along the x-axis, exciting the spin 
Hall torque pointing in the ^-direction on the magnetization 
m in the ferromagnet. The applied magnetic field is denoted 
as Ht. (b) Schematic view of the precession trajectory of 
the magnetization on the constant energy curve. The solid 
circle is the trajectory in the absence of the magnetic field 
or in the presence of the field along the 2 -axis, whereas the 
dashed elliptical lines are those in the presence of the field in 
the x and y- axes. The solid and dotted arrows represent the 
directions of the spin Hall torque and the damping torque, 
respectively. 

should be satisfied to excite the self-oscillation. If this 
condition is not satisfied, the magnetization directly 
moves to a constant energy curve including the saddle 
point without showing a stable steady precession, and 
stops dynamics because the spin torque does not balance 
the damping torque for E ln - m < E < E sa ddie■ An example 
of such dynamics is shown below; see Fig. [© We empha¬ 
size that Eqs. © and © are applicable to any kind of 
physical system showing a self-oscillation. 


III. SPIN HALL SYSTEM 

Let us apply the above discussions to the spin Hall sys¬ 
tem schematically shown in Fig. |T](a), where the electric 
current flows in the nonmagnet along the x direction, 
whereas the ferromagnet is attached along the 2 direc¬ 
tion. The spin polarization of the spin current is geomet¬ 
rically determined as p = e y . In the spin Hall system, 
the spin torque strength H s is given by 

Mj 

s 2 eMd' 


( 7 ) 
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where i9 and d are the spin Hall angle and the thickness 
of the ferromagnet, respectively. The magnetic field H 
consists of the applied field H t and the perpendicular 
anisotropy field HKm z e z . We can assume that H t > 0 
without losing generality because the sign of H t only af¬ 
fects the sign of j(E) derived below. Since we are in¬ 
terested in a perpendicular ferromagnet, we assume that 
Hk > H t > 0. Figure |T] (b) schematically shows the 
precession trajectory of the magnetization on a constant 
energy curve, where the directions of the spin Hall torque 
and the damping torque are represented by the solid and 
dotted arrows, respectively. The spin Hall torque is par¬ 
allel to the damping torque for m y > 0, whereas it is 
anti-parallel to the damping torque for m y < 0. This 
means that the spin Hall torque dissipates energy from 
the ferromagnet when m y > 0, and supplies the energy to 
the ferromagnet when m y < 0. Then, due to the symme¬ 
try of the trajectory, the net energy supplied by the spin 
Hall torque, W s , is zero when the applied magnetic field 
points to the x- or ^-direction. This means that Eq. © 
cannot be satisfied, and thus, self-oscillation cannot be 
excited in the spin Hall system in the absence of the ap¬ 
plied magnetic field, or in the presence of the field point¬ 
ing in the x- or ^-direction. Therefore, in the following 
we focus on the applied magnetic field pointing in the 
y-direction. The magnetic field and the energy density 
are given by 


m y (E) = r 3 + (r 2 - r 3 )sn 2 (it, k), (13) 


m z (E) = yjl - r§ - (r 2 - r|)sn 2 (u, k), (14) 


where u = qy 7 H t H\^/2^r\ — r 3 t , and rq are given by 

E 


n(E) = - 


MH t ’ 


(15) 


r2 (E) = — + Jl+(—'\ +2EL (16) 

V ' H k V MH k V ’ 


, x H t f H t \ 2 E 

Tz(E) ~ m ~ V 1 + \Th) + m7/k (I7) 

The modulus of Jacobi elliptic functions, sn(u, k) and 
cn(u, k), is 


k = 




r 2 - r 3 
r\ - r 3 ' 


(18) 


The derivations of Eqs. (El), El), and (1141) are shown in 
Appendix A. The precession period is 


H — Ht& y -f- Fk^®2) (^) 


E = —MHtrriy 


MH k , 
—^— m „ 


The minimum energy of Eq. © is 


E • — 


MH k 

2 



(9) 


( 10 ) 


which corr esponds to a point m s t a bie = 
(0, Ht/HK, -y/1 — (H t /HK) 2 )- On the other hand, 
Eq. © has a saddle point at m sa ddie = (0,1,0), 
corresponding to the energy density 

^saddle = —MH t . (11) 


Since the magnetization initially stays at the minimum 
energy state, and the magnetization dynamics stops 
when m reaches the saddle point m sa ddie) we consider 
the energy region of E min < E < E s ad die- To calculate 
Eqs. © and ©, it is necessary to solve a nonlinear equa¬ 
tion dm/dt = —qmxH, which determines the precession 
trajectory of m on the constant energy curve. Since the 
constant energy curve of Eq. © is symmetric with re¬ 
spect to the yz- plane, it is sufficient for the calculation of 
Eqs. © and © to derive the solutions of m for half of 
the trajectory in the region of m x > 0, which are exactly 
given by 


( 12 ) 


r(E) 


2K (k) 

1 ^/H t H K / 2 y/n - r 3 ’ 


(19) 


where K(fc) is the first kind of complete elliptic integral. 
The work done by spin torque and the dissipation due to 
damping, W s and #(*, are obtained by substituting Eqs. 
El): ED, and m into Eqs. © and ©, integrating 
over [0,r/2], and multiplying a numerical factor 2 be¬ 
cause Eqs. ED: ED: and ED are the solution of the 
precession trajectory for a half period. Then, and W a 
for E min < E < .Buddie are exactly given by 


_ 8 MH sy /n - r 3 

Fr § — - 

ZHtsjH^I(2H t ) 


( 20 ) 


= 

3y/H K /(2H t ) 

where Hs and H a are given by 


H s =H t 


1 — r 2 
ri - r 3 


K(A0- 


(~ 

\M 


HT) 

Hk) 


( 21 ) 


E (k), (22) 


(23) 

Here, E (k) is the second kind of complete elliptic inte¬ 
gral. The derivations of Eqs. ED and ED are shown in 


m x (E) = (r 2 - r 3 )sn(u, k)cn(u, k), 
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FIG. 2: The dependence of the current j(E), Eq. (124f) . 
for several values of Ht/H k on the energy density E. For 
simplicity, the horizontal and vertical axes are normalized as 
j{E)/j c and E / (E sad die F ln i T! ) [-£'min/(-£'saddle Trnin )] to 
make j(E m , n ) — 1, E m , n — 0, and E sad die — 1* 



FIG. 3: Typical magnetization dynamics excited by the 

spin Hall effect. The parameter values are taken from 
experiments ^ -38 ' 42 as M = 1500 emu/c.c., Hk = 540 Oe, 
a = 0.005, 7 = 1.764 x 10' rad/(Oe-s), d = 1 nm, § = 0.1, 
and H t = 50 Oe. The current magnitude is 14 x 10 6 A/cm 2 , 
while the critical current, Eq. (HU), is 13 x 10 6 A/cm 2 . 


Appendix B. The current j(E) for E m ; n < E < Caddie 
is given by 


2aeMd H t H a 

m m7’ 


(24) 


The currents for E E m [ n and E —> Saddle are^i 


^ _ 2 aeMd H K 

min)- M H J Hk 




(25) 


, 2aeMd {3H K -2H t \ 

= ~M~ { - 5 -)• (26) 

Equation (HHl) is the current density satisfying Eq. ©, 
or equivalently Eq. ©• Then, let us investigate whether 


Eq. (l24l) satisfies Eq. ©. It is mathematically difficult 
to calculate the derivative of Eq. (IM1) with respect to 
E for an arbitrary value of E, although we can confirm 
that j(E min ) > j{E sadd ie) for H t < H K . We note that 
a parameter determining whether Eq. © is satisfied is 
only Ht/ElK because the other parameters, such as a and 
M, are just common prefactors for any j(E). As shown 
in Fig. [2j j(E) is a monotonically decreasing function of 
E for a wide range of H t /H k, i.e., Eq. © is not satis¬ 
fied. This result indicates that the magnetization stays 
in the equilibrium state when j < j c = j(E m i n ), whereas 
it moves to the constant energy curve of E sadd i e without 
showing stable self-oscillation when j > j c because the 
spin Hall torque does not balance the damping torque 
on any constant energy curve between E m i n and E sa ddie- 
The magnetization finally stops its dynamics at ±m sa ddie 
because all torques become zero at these points. Figure 
[3] shows a typical example of such dynamics, in which the 
time evolution of each component is shown. Therefore, 
self-oscillation solely by the spin Hall torque cannot be 
excited in the perpendicular ferromagnet. This is a possi¬ 
ble reason why the self-oscillation has not been reported 
yet. 

Recently, many kinds of other torques pointing in 
different directions or having different angular depen¬ 
dencies, such as field-like and Rashba torques, have 
been propose d 28 ' 29 ' 36 ' 37 ' 40 ' 43- —. These effects might 
change the above conclusions. Adding an in-plane 
anisotropy 21,2 -, tilting the perpendicular anisotropy^, 
or using higher order anisotropy might be another 
candidate. Spin pumping is also an interesting phe¬ 
nomenon because it modifies the Gilbert damping 
constant^ - —. It was shown in Refsj 48 ' 50 that the en¬ 
hancement of the Gilbert damping constant in a fer¬ 
romagnetic /nonmagnetic/ferromagnetic trilayer system 
depends on the relative angle of the magnetization. This 
means that the Gilbert damping constant has an angular 
dependence. In a such case, it might be possible to sat¬ 
isfy Eqs. © and © by attaching another ferromagnet 
to the spin Hall system and by choosing an appropriate 
alignment of the magnetizations. The above formulas 
also apply to these studies. In Appendix C, we briefly 
discuss a technical difficulty to include the effect of the 
field-like torque or Rashba torque. 


IV. CONCLUSION 

In conclusion, we developed a method for the nonlinear 
analysis of the LLG equation in the spin Hall system with 
a perpendicular ferromagnet. We summarized physical 
conditions to excite self-oscillation by the spin transfer 
effect. The first condition, Eq. ©, or equivalently Eq. 
©, implies that the energy supplied by the spin torque 
during a precession on a constant energy curve should 
equal the dissipation due to damping. The second con¬ 
dition, Eq. ©, implies that the current to balance the 
spin torque and the damping torque in the self-oscillation 
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state should be larger than the critical current to desta¬ 
bilize the initial state. By solving the nonlinear LLG 
equation, we derived exact solutions of the energy sup¬ 
plied by the spin transfer effect and the dissipation due to 
damping, and showed that these conditions are not sat¬ 
isfied. These results indicate that self-oscillation cannot 
be excited solely by the spin Hall torque. 

The author would like to acknowledge T. Yorozu for 
his great constructive help on this work. The author also 
thanks M. Hayashi, H. Kubota, and A. Enrura for their 
kind supports. This work was supported by JSPS KAK- 
ENHI Grant-in-Aid for Young Scientists (B) 25790044. 


We note that Eqs. m, (ini). and m are periodic 
functions with the period given by Eq. (fTUl) . On the 
other hand, when E = Twaddle. the magnetization stops 
its dynamics finally at the saddle point m = (0,1, 0). The 
solution of the constant energy curve of E s addle with the 
initial condition m v { 0) = r 3 can be obtained by similar 
calculations, and are given by 


/ H t \ tanh(jzf) 
\ Hk ) cosh(z d) ’ 


(A6) 


m v = —1 + 


2 H t 
H k 



tanh 2 (j4), 


(A7) 


Appendix A: Precession trajectory on a constant 
energy curve 

Here, we show the derivation of Eqs. ED. m, and 
El- The precession trajectory on a constant energy 
curve is determined by dm/dt = -7m x H. The y- 
component of this equation is dm y /dt = ^HKm x Tn z . 
Thus, we find 


m z (l-—') -, (A8) 

y Hk V Hk J cosh (vt) 

where v = 7^H t (H K - H t ). 

Appendix B: Derivation of Eqs. (120ft and (121ft 


/ 


dt = 


1 

jHk 


duly 

m x m z ' 


(Al) 


As mentioned in Sec. mu since the constant energy curve 
of Eq. © is symmetric with respect to the yz- plane, it 
is sufficient to derive the solutions of m for half of the 
trajectory in the region of m x > 0. Using E and m v , m x 
and m z are expressed as 


m x 


1 - m 2 + 


2 E 
MH k 


2 H t 
H k 


?7l 


y> 


(A2) 


m z 


2 E 


2H, 


MH k H k 


Oy. 


(A3) 


The initial state of m y is chosen as m y ( 0) = r3, where 
r 3 is given by Eq. El- Then, m y at a certain time t is 
determined from Eq. (ED as 


1y/2H t H K / dt 
Jo 


dm ' 


(A4) 


\J(™' y ~ n)(m' y - r 2 )(m' y - r 3 ) 


We introduce a new parameter s as m y = ?'3 + (r 2 — ^3)s 2 - 
Then, we find 


= J 


ds' 


0 \J{1 - s ,2 )(l - k 2 s' 2 ) 


(A5) 


where the modulus k is given by Eq. m- The solution 
of s is s = sn(u, k). Therefore, m v is given by Eq. (fl3l) . 
Equations El and (flTl) are obtained by substituting Eq. 
El into Eqs. (IA2I) and (IA3I) . 


Using Eqs. El. El. and El. the explicit form 
of Eq. I© for the spin Hall system is given by W s = 

7 M H s f dtw s , where w a is given by 

w s = (H t - H K r 3 ) (1 - r\) 

+ {-2 H t r 3 + H k [r 3 (r 2 + r 3 ) - (1 - r 2 )] } (r 2 - r 3 )sn 2 (u, k) 
+ {-H t + H K (r 2 + r 3 )} (r 2 - r 3 ) 2 sn 4 (u, k). 

(Bl) 

Similarly, Eq. m for the spin Hall system is given by 
= — cryM / dtw a , where w a is given by 

w a = ( 1 - r\)(H t - H K r 3 ) 2 

- [2 H 2 r 3 - Hl(r 2 + r 3 )(l - 2r 2 ) + 2H t H K (l - r 2 r 3 - 2r 2 )] 

x (r 2 — r 3 )sn 2 (u, k ) 

- [H t - H K (r 2 + r 3 )f (r 2 - r 3 ) 2 sn 4 (u, k). 

(B2) 


Then, W s and W a are obtained by integrating over 
[0,t/2], and multiplying a numerical factor 2. The fol¬ 
lowing integral formulas are useful, 



u — E[am(u, k), k] 
k 2 


(B3) 



sn(u, fc)cn(u, k) dn(u, k) 
'' 3 k 2 

2 + k 2 
+ ^ U 
2(1 + k 2 ) 

3 fc 4 


E[am(u, fc), fc], 


(B4) 


where E(u, fc), am(rt, fc), and dn(tt, fc) are the second kind 
of incomplete elliptic integral, Jacobi amplitude function, 
and Jacobi elliptic function, respectively. 
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Appendix C: The effect of the field-like torque or 
Rashba torque 

The direction of the field-like torque or the Rashba 
torque is given by m x p, where p is the direction of the 
spin polarization. This means that the effects of these 
torques can be regarded as a normalization of the field 
torque m x H. Then, the energy density E and the 
magnetic field H in the calculations of W s and W a should 
be replaced with an effective energy density £ and an 
effective field 23 given by 

£ = E - /3MH S m ■ p, (Cl) 


B = H + (3H s p , (C2) 

where a dimensionless parameter /3 characterizes the 
ratio of the field-like torque or Rashba torque to the 
spin Hall torque. We neglect higher order terms of the 
torque^ for simplicity, because these do not change the 
main discussion here. In principle, j(E) satisfying Eq. 
© can be obtained by a similar calculation shown in 
Sec. mo However, for example, the right-hand-side of 
Eq. (l24l) now depends on the current through £ and H. 
Thus, Eq. (IMll should be solved self-consistently with 
respect to the current j, which is technically difficult. 
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